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The cross section for 77 — > 7r°7r° and the pion polarizabilities are computed, within 
generalized chiral perturbation theory, in the full one loop approximation, i.e. up to and 
including order 0(p 5 ). The result depends on the parameter ot^ defining the tree level ir — 71 
scattering amplitude and on an additional low energy constant. The latter is shown to be 
related by an exact sum-rule to the e + e~ data. The parameter a n7T is related to the quark 
mass ratio r = m s /rh via the expansion of pseudoscalar meson masses. The generalized one 
loop 77 — > tt°7t° amplitude agrees with the experimental data in the threshold region provided 
r = m s /m < 10. Higher order corrections are estimated comparing our calculation with the 
dispersive approach. 
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1 Introduction 



The reaction 77 — > 7r 7r at low energies provides a so far rare example of a probe of 
the chiral symmetry breaking sector of QCD for which experimental data already exist H. 
Moreover, these data do not exhibit any unexpected feature, as they can be reproduced with- 
out difficulty using the standard methods based on dispersion relations, unitarity and 
resonance saturation. The one-loop prediction of the standard chiral perturbation theory M 



(xPT) appeared prior |], [K| to the publication of the experimental results. It does not involve 
any free parameter and it disagrees with the data by several standard deviations even close to 
threshold. Recently, the full two-loop calculation has been completed by Bellucci, Gasser and 



Sainio [11]. It involves three new 0(p 6 ) constants (estimated via resonance saturation) and 



it agrees with the experimental data within errors. The purpose of this work is to add one 
more information to this list: In the framework of generalized x?T (defined in ]T2| , |TB|), the 



agreement with experiment near threshold is already reached within the one-loop approxima- 
tion, provided the ratio of quark masses r = m s /m is considerably lower than usually expected, 
typically r < 10. 

The important question to ask is: What can 77 — ► 7r°7r° data tell us about the expansion 
of the QCD effective Lagrangian [|14|, 15| ? While each individual term of C e ^' is uniquely 



determined |16| by symmetry properties of QCD, the relative importance of different terms 
depends on the actual values of the quark masses and of the low-energy constants characterizing 
the chiral structure of the QCD vacuum. The latter in turn determines how the expansion of 
£ e ff- in powers of quark masses and external momenta should be organized such as to obtain 
a good convergence rate. The standard chiral perturbation theory || is an expansion which 
assumes that the quark masses m u , and m s are small enough not only with respect to the 
hadronic scale Ah ~ 1 GeV, but also as compared to the scale Bq of the single flavour quark- 
antiquark condensate, < qq >o= —BqFq, defined at m u = ma = m s = 0. Within QCD, the last 
assumption is hard to justify a priori, since nothing prevents Bq to be as small as, say, the pion 
decay constant F ~ 90 MeV, i.e. much smaller than A#. To illustrate the importance of the 
condition m q <C B , let us mention the fact that the standard relation between pseudoscalar 



meson and quark mass ratios ]17 



m M 2 1 

-^ = 2— f -l + --- = 25.9 + --- > rh = -(m u + m d ) , (1) 
m Mi 2 



receives, at the next order, corrections || (represented by the ellipses) which are of both 
types, 0{m q / Ajj) and O(m q /B ). On the other hand, the analysis of the Dashen-Weinstein 
sum-rule for the deviations from the Goldberger-Treiman relation suggests that r might differ 
considerably from the above value [19|, e.g. r < 10. If confirmed^, this result could be in- 



terpreted as an experimental indication of an important 0(m q / Bq) contribution to Eq. ([I]). 

1 For a recent determination of the 7rN coupling constant, see Ref. [5| 
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The generalized chiral perturbation theory [|12], [13| reformulates the expansion of the effective 
Lagrangian without assuming that m 9 -C _B . It is as systematic and consistent as the standard 
%PT. At all orders in p/ Ah and in m 9 /A#, both expansion schemes, standard and generalized, 
sum up the same effective Lagrangian C e ^' of QCD. However, at each finite order, the gener- 
alized xPT takes into account terms which standard relegates to higher orders. Already 
at the leading 0(p 2 ) order, the generalized involves more terms and more parameters: In 
particular, the quark mass ratio r is a free parameter, to be determined by experiment. In the 
limit r = r 2 = 2M\jM 2 — 1, one recovers the standard as a special case. In practice, the 



generalized is defined by the formal counting rule rh, m s , B ~ 0(p) [|1^, [13]], replacing 
the usual rule rh, m s ~ 0{p 2 / Ah), B ~ 0(A#) characteristic of standard [[|. Hence, 
for p <C Ah, one can write 

C eff. = C (2) + £ (4) + £ (6) + . . . = £(2) + £(3) + £ ( 4) + £(5) + £(6) + . . . ? (2) 



where the first expansion is the standard one |J, whereas the second one corresponds to the 
generalized x?T. 

In the next section, we summarize our results for the one-loop generalized 77 — > 7r 7r 
amplitude. It depends on the generalized tree-level tc—tt amplitude |T2|, |TB[ (kaon loops give only 
a small contribution) and it receives a constant shift from an order 0(p 5 ) tree- level contribution. 
The latter is shown to be calculable, through a low-energy theorem, from e + e~ — *> hadrons data 
(Section 3). As an indicative estimate of higher order corrections, we submit, in Sect. 4, our 
result to a dispersive analysis along the lines of Ref. 0. We finally add a few comments on 
charged and neutral pion polarizabilities. 

We find it convenient to perform the whole analysis within the three-light-flavours %PT. 
It allows us to keep kaon loops under control and, mainly, it should help in establishing re- 
lationships to other observables (c.f. Sect. 3) and to similar processes such as rj — > 7r°77 



2 77 —> 7T 7T to generalized one loop order 

In the present section, we establish the expression at generalized one loop order for the 
amplitude of the reaction 77 — > 7r°7r°. The modified chiral counting specific to generalized 
%PT leads to the appearance of terms of odd orders, which do not correspond to an increase 
in the number of loops, but to additional corrections in powers of the quark masses. Thus, 
the leading order contributions from C^ 2 ' receive tree- level corrections from C^ 3 \ which come 
before the one loop effects. Similarly, before going to the two loop order, one has to take into 
account contributions at order 0(p 5 ). The generating functional (S-matrix) at one loop order 
in generalized is thus given by the following expansion: 

7 eff. _ 5(2) , 5(3) , 5(4) , 5(4) , 5(5) , 5(5) , (o\ 
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Here, Z\ ree corresponds to all the tree diagrams made from an arbitrary number of vertices from 
£( 2 ); Z^le denotes the same set of diagrams, but where a single vertex from is replaced by 
a vertex from C^ 3 \ Z^\ e contains a single vertex from or two vertices from , with an 
arbitrary number of additional vertices from £( 2 \ Similarly, Z^\ e stands for all tree diagrams 
made from any number of vertices from and either a single vertex from or one vertex 
from and one vertex from . Finally, z[fl op stands for all one loop diagrams made with 
vertices from only, while Z^ oop denotes the same diagrams, but with one of those vertices 
replaced by a vertex from 

The general analysis of Z e ^- up to and including the order 0(p 5 ) (Eq. (|)) and of the 
various low-energy constants involved is beyond the scope of the present work, and will be given 
elsewhere [^2|0 . Here, we shall merely consider those terms which contribute to the matrix 
element 

< 7i ( Pl )7r°(p 2 ) out\ 7 (A; 1) ei ) 7 (fc 2 , e 2 ) in >= ze 2 (2vr) 4 5 4 (P / - P^M^ , (4) 

withQ 

= H(s, t, u){-n^ - k lu k 2f i) H • (5) 

(The ellipses stand for the other Lorentz tensors in the general decomposition of M^, whose 
form factors receive contributions, for on-shell photons, only from order 0(p & ) onward.) In 
this case, the terms z£l e , Ztrle an d Ztrle m Eq. ([|) do not contribute. The remaining three 
terms contribute to the amplitude H(s, t, u), the loop diagrams z[fl op and Z^ and the tree 
diagrams of order 0(p 5 ) are all separately finite. Topologically, the loop diagrams are the same 
as in the standard case [[|, ||10|| , but the vertices are to be read off from and from C^ 3 \ 
where 

£(2) = ?ll( D »u + D fl U}+2Bo{M q {U + C/+)) 

+ A Q {{M q U) 2 + (M q U + f) (6) 
+ Zg(M q U + M q U + ) 2 + ..-} , 

£( 3 ) = ^UiD^U+D^UiMqU + U+M q )) 

+ ^(D^U+D^U) (M q (U + + U)) 

+ Pl ((M q U) 3 + (M q U + f) + P2 (M 3 q (U + U + )) (7) 

+ P ±((M q U) 2 + (M q U + ) 2 )(M q (U + U+)) 

+ P5 (M 2 q )(M q (U + U + )) + P7 (M q (U + U + )) 3 + -- 



and 



2 For a discussion of the leading order, see Ref. [Q 

3 The amplitude denoted by H(s, t, u) is called A(s, t, u) by the authors of Ref. [JlTJ . We have adopted a 
different notation in order to avoid confusion with the standard notation A(s 1 1, u) for the ir — ir scattering 
amplitude. 
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with at most one vertex per diagram coming from £( 3 \ The notations are as given in Ref.||, 
Aiq = diag(m, m, m s ) denotes the quark mass matrix and the ellipses stand for those terms 
which do not contribute to H(s, t, u). The tree graph contributions of order 0(p 5 ) are described 
byQ 



£(5) = c 1 {F^M q F Rfll/ U + F R iV M q F LfMU U + ) 

+ c 2 (FrU + F RlMU UM q U + F^UF Lflu U + M q U + ) 

+ c 3 (F^F R , u (UM q + M q U + ) + F£ U F L ^(U + M q + M q U)) 

+ c 4 (Fjf + F^F L ^)(M q {U + U + )) + ... . 



(8) 



The final result will only depend on F n , M n , and on a few other independent combinations of 
rh, of r, and of the various low-energy constants occurring in (|6|), (|7|) and (|8|). The relevant 
vertices are given in Fig. 1. Up to the order 0(p 5 ), the amplitude H(s, t, u) depends only on 
s, H(s, t, u) = H(s) + 0(p 6 ), and takes the following form 



H(s) = A^^is) + A K - loop (s) - i£ 

Fz 



(9) 



The pion loop contribution A n loop may be written as 

4 



.4 



TV — loop t 



^M 2 

3 



Fl 



+ a„ 



^\Gi^ 
3Ff J K Ml 



(10) 



where G denotes the following loop function 

l + ~( ln I^+^J 
l-farctg 2 (^ 



16w 2 G(z) 



if 4 < z 
if < z < 4 
if ^ < 



= Vl - 4/* 



(11) 



Here and are two combinations of the low-energy constants of CS 1 ' 
parametrize the 0(p 3 ) tree-level tc — tt scattering amplitude |T2| 



A(s \t,u)= p m 3 2 *- + a n7T ^ + 0{p A ) . 



which 



(12) 



The kaon loop contribution has a similar structure, 



A K-loap, s s 



[(M* - il^) 2 + IMvMKOvjh 



Ml" 



(13) 



4 In standard xPT, these terms would be part of the 0(p 6 ) contributions ITTf 
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in terms of two other combinations of the low-energy constants, a w K and fi^K-, which parametrize 
the 0(p 3 ) tree-level K — n amplitude A + (s, t, u) f23|, symmetric in s and u: 

*> u ) = §§ (* - \ M l ~ \ M k) + ^2 - M -) 2 + 2M.M^«^] + 0(p 4 ) . (14) 
It is useful to express a W7r and ct^ as 

a wK = a l ^ d + da^K , (15) 



where a l ^ d and a l ^jf arise from C^ 2 \ whereas 5a nn and 5a n K stem from £( 3 ). The explicit 
expressions for 5a ww and for 5a w K in terms of the low-energy constants appearing in and 
in C^- 3 ' will not be used in the sequel and can be obtained from the formulae collected in the 
Appendix. Similarly, one finds: 

A™ = 1 + 2£m + 4|m , = 1 + 2£m + 4£m(3 + r) . (16) 

(In both cases, the leading order values are f3 l ^ d = = 1.) Finally, we come to the tree 
contribution from which gives a constant shift to H(s), see (|9|), with 

10 / n w 16 . 10, n _ 

c = — — (ci + c 2 + 2c 3 )m - yc 4 m « — — (ci + c 2 + 2c 3 )m , (17) 

where the second, approximate, expression is obtained by invoking the Zweig rule. This constant 
shift in H(s) appears (along with other contributions) only at order 0(p e ) in the standard 
approach and c is in fact related to the constant d 3 of Ref. fi~H| ( see the discussion at the 
end of Section 3). However, since the loop contributions to H(s) are finite by themselves, c is 
not renormalized at the order we are considering The value of c depends on r and will be 
discussed below. 

First, let us consider the amplitude H(s) when contributions up to order 0(p 4 ) alone 
are taken into account. In that case, the constants c, 5a nn and Sa-^x drop out, and only the 
leading, 0(p 2 ), expressions for the n — n and K — tt amplitudes are involved. These essentially 
depend on F n , M n and r, since Hl2| , |D| , p3| , 



^ = 1 + 6-1^(1 + 20 , Pl e : d = l, r 2 = 2^f -1-25.9, 

a^ = l + ——(a l *f-l) , # = 1, r x =. 2-/-1-6.3, (18) 
ri + 1 M n 



5 In terms of the K — tt isospin amplitudes, A + — |A 3 / 2 + \A X ^. 

6r 



°The fact that c is finite in our case, whereas the constant of is renormalized does not bear any 
contradiction, and will be explained later (Sect. 3). 
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and since the Zweig rule violating parameter £ = /A is expected to be small. In fact, 



vacuum stability arguments impose the following bounds at leading order [JL3 

r > 7*1 , 



(19) 



and 



o<C< 



1 r — 7*i r + r x + 2 

2 r 2 — r r + 2 



(20) 



a 



lead 



„lead 



. For the critical 
4 fl2fl, |E3|, and for an intermediate value, say r ~10, 



For r = 7*2, one recovers the standard values of Weinberg [[17 
value r = 7*1, one obtains = ®nK 
2 < a^" d < 2.3. Thus, as r decreases from its standard value r 2 towards its critical value r 1; 
u l ™ d varies by a factor four. This results in an increase of at order 0(p A ) in the threshold 

region of about 25% when r ranges through the same values. The corresponding variation of 
the cross section is shown in Fig. 2. 

At full one loop order, i.e. taking into account the order 0(p 5 ) corrections, the amplitude 
H(s) no longer depends on r alone. It also involves the low-energy constants of £® and the 
combination (|i~7D of low-energy constants from . We first notice that the constant £ involved 
in the expressions (|T6|) for (3^ and (3^k is also responsible for the splitting between F n and Fx 
at order 0(p 3 ) [jllj (in the standard case, this splitting occurs only at order 0(p 4 )): 



Fl = F 2 [1 + 2m£ + 2ml{2 + r) 



(21) 



F 2 K = F 2 [1 + m£(l + r) + 2m£(2 + r) + 



Therefore, 



1 



F 2 

% - 1 I + 



r - 1 V Fl 

TV 

which allows us to express (3 nn and fl^K in terms of r 

2 



1 + 
1 + 



Ft 

r-1 VF2 

-l U 2 



r 



-1 (1 + 2^/0 



1 [1 + 2(3 + r)£/£] 



(22) 



(23) 



(24) 



up to the ratio which is anyhow expected to be small, due to the Zweig rule. Notice that £ 
and £ are related to the low-energy constants L 5 and L 4 of Ref . [Q , respectively. Up to higher 
order, 0(p 6 ), corrections we may further replace r in the above formulae by its expression 



r(a nn ) obtained from Eq. (|l8|) , but with a l ™ d replaced by = 
= 1.22). Finally, we establish, in Sect. 3, the following sum-rule: 

5 l r r°° ds 



a 



lead 



+ 5a W7r (we take F K /F n 



1447T 5 



— \R + (s) - 3RJs)] --\n(^ 

AMI S 1 +K ' V n 2 



+ 0(m 2 m q B, 



(25) 



6 



where R + (R-) represents the cross section for e + e~ into hadronic final states with even (odd) 
total G-parity, normalized to the cross section for e + e~ — > Again, r in the formula above 

is to be thought as expressed in terms of a n7T . The evaluation of the integral in the narrow 
resonance approximation gives the following estimate: 

' (4.6 ±2.3) x 1(T 3 . (26) 



r — 1 



Hence, up to the small contribution from the kaon loops, and up to the Zweig rule violations 
contained in ( and in we have been able to express the cross section for 77 — > 7r°7r° at 
order 0(p 5 ) in terms of the single parameter a n7T . This constant could, in principle, be extracted 
from ix — ix phase shifts, given sufficiently accurate datafj. We recall that if standard is 
correct, then a n7r is close to 1. The smaller B , the more a nn increases towards a W7T ~ 4. 
In Fig. 3, we show the cross section for the three values a nn — 1, 3 and 4. For a n7T > 3, 
our result agrees, within errors, with the data in the threshold region. For a nn = 4 (i.e. 
when the factor l/(r — 1) in and ( P5| ) is largest), we have checked that the cross section 
is indeed insensitive the the contribution coming from the kaon loops, by varying a n K in the 
range 2 < a n K < 6, and by taking /3 n K from Eq. (|24]) . We have also checked that the result is 
insensitive to the variation of the constant c within the range —1.30 x 10~ 3 < c < —0.44 x 10 -3 
given by Eq. (|26|) for a nn = 4. Actually, even if we set c equal to zero, the low-energy cross 
section is only barely affected, see Fig. 4. On the other hand, the cross section is more sensitive 
to variations of /3 nn coming from different values of the Zweig rule violating ratio £/£. Taking 
this ratio in the range —0.2 < £/£ < +0.2, which amounts to a Zweig-rule violation of 40% in 
(and corresponds to the range allowed for L 4 /L 5 in [|J), we find the variations in the cross 
section as shown in Fig. 5. 



3 A low-energy theorem for the 0(p 5 ) constant c 

In this section, we derive the sum-rule ( p5|) and discuss its numerical evaluation. Consider 
the decomposition of the electromagnetic current into its 1=1 and 1=0 components: 

,3 1 



J, = V£ + , (27) 

where 

K = \{ui»u - d^d) , = ^(«7m m + dj^d - 2s^^s) . (28) 

The corresponding vacuum polarization functions, 

i J d 4 xe ! ^' < n\T{V«(x)V u a (0)}\n >= [q^q u - ^ u q 2 )W a (q 2 ) , a = 3, 8 , (29) 
7 An analysis [Fl3| of the available data shows that any value of in a range ^1 — 4 can be accounted for. 
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and 

i J d'xe^ < n\T{j;(x)f u (o)}\n >= fog* - v»vq 2 W- m iq 2 ) , (30) 



satisfy 



n e - m -(Q 2 ) = n 3 V) + ^V) , (31) 



if isospin violations due to the electromagnetic interactions and to m u ^ are neglected. 
Furthermore, the spectral density of the electromagnetic current is related to the inclusive 
e + e~ cross section 



-^mW- m (q 2 ) = ~^R(q 2 ) , (32) 

7T 1Z7T Z 



where, as usual, 



R(q 2 ) = °^^f] = R + (q 2 ) + R_(q 2 ) , (33) 



and where, corresponding to the decomposition (p7|), R(q 2 ) has been split into its contributions 
from G-parity even and from G-parity odd hadronic final states, with 

-QmU 3S (q 2 ) = ^- 2 R + {q 2 ) , -3mn 8 V) = -^R-tf) ■ (34) 

The difference IT 33 (g 2 ) — II 88 (g 2 ) then satisfies an unsubtracted dispersion relation, and conse- 
quently one can write the following sum-rule: 

n 33 (o) - n 88 (o) = -f - [3mn 33 ( s ) - s m n 88 ( s )] . (35) 

7CJ 4M2 S 

On the other hand, the left-hand side of Eq. (^) above satisfies a low-energy theorem: From 
the generating functional Z e ^' at order 0(p 5 ), one obtains: 

n 33 ( S ) = - 2[L{M + 2H r M) 

- 4(ci + c 2 + 2c 3 )m - 8c 4 m(2 + r) 



l (s -4MM4) + ^-4<)J(^) (36) 



2 1 /, Ml \ 1 1 / M 2 K \ 1 „ 



3 32tt 2 V fi 2 J 3 32tt 2 V /i 2 / 48vr 

n 88 (s) = 2[2^(Ai) + 2flT(Ai)] 
4 

- -(1 + 2r)(ci + c 2 + 2c 3 )m - 8c 4 m(2 + r) 

+ h s - iM «^ (37 > 



.V 2 - \ I 



32tt 2 V h 2 J 48tt 



In + 1 + -— + 0(s, m 2 , mB, 



The low-energy constants L r 1Q (n) and H[(n) are as defined in ||: they are common to both 
£( 4 ) and since they are not associated with symmetry breaking terms, 

L 10 (U + F Rflu UFn + HxiFB^Fj? + F Lflv F£ v ) e £ (4) , £ (4) . (38) 

They depend on the renormalization scale /x, again as specified in [g], but both II 33 and n 88 
are independent of fi. The loop function J reads 

16tt 2 1 v l + a 



J ( z ) = 7ZZ*{*( l *TTZ: + i «) +n » <r(z) = y/l-*/z, (39) 



for z > 4, and has the following expansion near 3 = 0, 

Taking the difference of Eqs. (f53) and (|57D , and neglecting the Zweig rule violating constant 
C4 as in Eq. (0), one obtains the following sum- rule: 

c ~ -y(ci + c 2 + 2c 3 )m (41) 



12 r - 1 



1 r ^ [3mn 33 ( S ) - Smn 88 ( S )] - -L In ( ^f] 
J 4M3 s K ' v ;J 24vr 2 UJ 



7T 



+ 0(rh 2 , rhBi 



which can now easily be written as given in Eq. Notice also the presence of the overall 

factor l/(r — 1) which suppresses the value of c when ~ 1. The above sum-rule is similar 
to and should be compared with the sum-rule relating the scale invariant combination 

L 10 = L[M + ^(lnM^/i 2 + 1) (42) 



to the difference of the spectral densities of the 1 = 1 vector and axial currents [pq] . 

The right-hand side of Eq. ( |4l| ) expresses c as a difference of two contributions. The first 
one, represented by the dispersive integral, is numerically the most important one and positive, 



so that c comes out with a negative value. Writing this contribution as in Eq. (|25j ) suggests 
a direct evaluation of the dispersive integral in terms of e + e~ data. However, such data are 
not so easy to obtain, especially in the isoscalar channel, and are affected by large error bars. 
We refer to Ref. [2f| for an overview of the experimental situation and for references. For our 



part, we shall evaluate the dispersive integral in the narrow resonance limit, taking into account 
contributions from p(768), u>(782) and 0(1020): 

1 ^m»( s| -«»(.]|„. = » {£^-3.5^-3 r ' 



7T J 4A/^ s w v n Ana 2 I M p M w M 

~ (11.1 ±2.0) x 10" 3 . (43) 

The experimental values for the e + e~ widths have been taken from the Particle Data Group 
compilation [[27j , and the error comes from the uncertainties on these numbers. The second term 
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on the right-hand side of Eq. (JO]) comes from contributions of virtual 2tc and IK intermediate 
states to the vacuum polarization functions II 33 (0) and II 88 (0). Numerically, we have 

1 i n ( Mf) = 5.4 x 1(T 3 . (44) 



24tt 2 V M 7 . 

We estimate the contributions of non-resonant states to the dispersive integral by the size of 
this chiral logarithm. Indeed, if we evaluate the dispersive integral ( [4"3"D replacing 
3mII 33 (s) — Smll 88 (s) by the imaginary parts as given by Eqs. fl36|) and fl3T|), we precisely 
obtain this logarithm: 



\ /I T |3mnS « " WW = ^ m (£) . (45) 



Hence, we do not add it to the contribution fl43|) , but take it as our estimate of the uncertainty 
on the value of c, which we thus obtain as 

c = - — j- ■ (4.6 ± 2.3) x 1(T 3 . (46) 

For a n7T = 4, one thus obtains c = —(0.87 ± 0.43) x 10~ 3 , whereas for a n7T = 1 this value 
comes out about five times smaller, c = — (0.18 ± 0.09) x 10~ 3 . This last value is relevant for 
the comparision with the constant d% of Ref. 11]. At first sight, the fact that the latter is 
renormalized, whereas our 0(p 5 ) constant is finite, might appear as puzzling. The point is that 
the divergent part of d 3 comes with an additional power of B Q , and thus counts as order 0(p 6 ) 
in generalized x?T. Therefore, the relationship between these two constants may be expressed 
as follows: 

9F 2 c 

d r M = -J^- q I 1 + O(B ln/i, m) + • ■ • } , (47) 

where the ellipses stand for additional corrections which might appear due to the fact that the 
constant c is defined in the three-light-flavours expansion scheme. In order to compare with the 
value of (ig extracted from standard x?T, in the numerical evaluation of Eq. ([IT]) above one 
should take the value of c for a n7T ~ 1, corresponding to the standard case, and consistently 
replace, up to higher orders, 2rhB by M%, i.e. consider B ~ 0(A#). This gives (notice that 
the sign is well defined) = (9.4 ±4.7) x 10 -6 + O(B ln/i) H — •, which is not incompatible 

with, say, the value of the scale invariant constant as computed from Appendix D of Ref. 



11 



d 3 = ^ n " s s ~ ±3.9 x 10" 6 . (48) 
3 8M§ 



It is possible to discuss the relationships between, say, £ and L 5 , or between £ and L 4 , along 
similar lines. This would, however, lead us too far astray, and we defer this discussion to 
forthcoming publications p8| , p2 . 
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4 Comparison with the dispersive approach 



Morgan and Pennington || have devised a dispersion relation treatment of the 77 — ► 7r°7r° 
amplitude which proved capable of representing the data quite correctly. The drawback of this 
type of approach is that the general principles of unitarity and analyticity do not completely 
constrain the amplitude. Donoghue and Holstein || argued that additional constraints should 
be provided by the chiral expansion. In this section, we discuss this idea, and we use the 
dispersive method, in conjunction with our 0(p 5 ) results, as a means of estimating the size of 
higher order chiral corrections^]. 

In order to implement the method, one has to consider the charged channel 77 — > 7t + tt~ 
as well as the neutral one. Let us designate by M xx , and M xx , the corresponding helicity 
amplitudes, normalized as in Eq. (3), and introduce the J=0 projections: 

fci*) = ^J d ^ M + + . /*(*) = ^J <^ M + + . (49) 

where /jv(s) is related to the amplitude defined in Eqs. (4), (5), H(s,t,u) = H(s) + 0(p 6 ), 
by H(s) = 4f N (s)/s. Next, one introduces isospin combinations fi(s), I — 0, 2 [|| 0] which 
diagonalize the elastic unitarity relation: 

fc = \h + \f2, /jv = |(/o-/ 2 ) ■ (50) 

Below the KK threshold, the following representation holds for these functions || : 

fj(s ) = l I (s ) + n T (s )s\ 7l + (s -s)- / — , nJ V V J \ \ , (51) 

I 7T JAM* x (x — s ){x — s — ie) J 



where is the Omnes-Mushkelishvili |29j function constructed from the n — n phase shifts 

Si(s) in the isospin channel /, 



fli(s) = exp 



s f°° dx 8i(x) 



71 J AMI x (x — s — ie) 



(52) 



(normalized such that fi/(0) = 1) and So is an arbitrary subtraction point. The function h{s) 
must have the same left-hand cut structure as the amplitude fi(s) and be analytic in the rest of 
the complex s plane. The representation ( |5TD then garantees that fi(s) has the correct analytic 
structure. In particular, the discontinuity along the right-hand cut 4M% < s < 00 is given by: 

Qmfi(s) = //(s)e-^ (s) sin <Jj(s) . (53) 

In the left-hand cut function one can separate the single pion pole (Born term) 

AM 2 s 
= ^( g /Ma arCtanhff( A^ ) + • (54) 



3 We are indebted to J. Gasser for repeatedly emphasizing this point to us. 
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The remaining piece, li(s), is generated by two or more pion exchanges in the t (u) channels 
and shows up in the chiral expansion starting at order 0(p 6 ). A factor of s has been separated 
in Eqs. (|5T| ) and (|54|) on account of low-energy electromagnetic theorems and a single arbitrary 
constant (in each isospin channel) appears due to boundedness properties for asymptotic values 
of s. 

It should be clear that a rigorous evaluation of the left-hand cut function li(s) is no 
simple matter. In fact, this can be done only order by order in chiral perturbation theory. For 
the purpose of doing order of magnitude estimates, however, it seems reasonable to retain the 



contributions of a few low-lying resonances: p, to, a\ and b\ [^OJ. The low energy region of 
interest here, 2M n < E nn < 450 MeV proves to be fairly insensitive to the details of which 
resonances are included and how they are parametrized |§. This leaves one with the task of 
determining the arbitrary constants 77. Following the suggestion of ||, one should be able to 
do so by matching the dispersive representation (|5l|) with the prediction, f?(s), in a point 
where the latter can be trusted. It is of course convenient to choose the subtraction point in 
( |5T| ) equal to the matching point, so that the constants 77 are simply given by: 

fi(so)-li(s ) (55) 

One has then to extend the calculation to the charged channel and one finds, for the 

isospin functions f*(s): 

™ = ^^'^) + ^w c ^) + ^w e ^ )+?i5 ^-^ (56) 

where 

A o = 7^ [Ml(5a nw - + , A 2 = 7^ Wl (2a™ + 4/U) - , (57) 

bri brf: 



A K - — — 

~ 4F 2 



p vK (~s - Ml - M 2 K ) + 2a 7rK M n M K + (M K - M n ) 2 



^2=0, (58) 



and where the function G was defined in Eq. (11). ff contains a contribution of order 0(p 2 ) + 
0(p 3 ) + 0(p 4 ) as well as a a set of contributions of chiral order five which are collected in 
the constants Aj. The combination c = 2(Ao — A2)/3 was evaluated in Sec. 3. In contrast, 
we do not know of a simple way to estimate the combination c + = (2A + A x )/3 which 
contributes to the tc + tt~ channel. In principle, of course, Aj can be evaluated from experiment 
via ;\TT analysis of the reactions 7r — > ewy, K — > ez/7 and of the mean charge radii of the pion 
and the of kaon. One expects the size of c + to be of order 10-30% that of L 9 + L w , so that it can 
be absorbed to some extent into the uncertainty in the numerical value of Lg + L%q (note that 
for consistency this numerical value itself should be determined here from an 0(p 5 ) analysis). 
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We have included the contribution of the kaon loop in (|5q) for completeness, but since the 
corresponding KK discontinuity is not accounted for in the dispersive formula (this could in 
principle be done following Ref. []3T]), we do not include it when we perform the matching. It 
is worth noting that the chiral formula (^) satisfies the general dispersive representation: at 
this order, one has to set h(s) = 0, tti(s) = 1 and 3m(l/0/(x)) = — a(s/M^)Af(s)/16TT in Eq. 



One should be aware that there are some ambiguities in the dispersive approach. Consider 
first the question of the matching point So in (fjo]). The authors of Ref. ||, for instance, 
have matched the 0(p A ) ;^PT amplitude with a dispersive representation not including the 
vector meson contributions at s = 0. A priori, one would expect that any point in the region 
< s Q < 4M^, where the amplitude is real, should be as good as another. In practice though, 
the results are sensitive to which point is chosen. We will show the results corresponding 
to the two extreme cases sq = and sq = 4M%. A further source of uncertainty, in the 
dispersive approach, comes from the n — tt phase-shifts themselves. Clearly, what will matter 
most in the calculation is how the phase-shifts behave close to threshold. This threshold 
behaviour is usually parametrized in terms of the scattering lengths and of the slope parameters. 
Experimentally though, the phase-shifts have been extracted only for energies larger than 600 
MeV. Extrapolation down to the threshold, even upon using constraints from the Roy equations 
and from KI4 data, is subject to some uncertainty. The influence of this uncertainty on the 
77 — > ir°7T° dispersive amplitude has been considered in ref. |§. Furthermore, in our case, the 
matching requires that one knows the relation between the constants a n7T , (3^ and the set of 
phase-shifts. This relationship has been studied in Ref. |13| and we will rely on its results in 



what follows. Consider the simple parametrization proposed by Schenk [j32|j : 

/ 4 /if 2 r - t s 



7T 1 



(59) 



where hi is related to the slope parameter bi by bi — bi — aj4M^/(E] — 4M^) + (a/) 3 . For 
I = we take two sets of parameters resulting from two fits to the production data obtained 
by fixing the value of the scattering length to a® = 0.26: 

a) a = 0.26, b = 0.203, c = -0.0126, E = 813.3MeV (60) 

b) a = 0.26, b = 0.324, c = -0.0274, E = 863.1MeV , (61) 



which correspond to the data of Eastabrooks-Martin [33] and of Ochs [B3J], respectively. The 
two sets differ essentially by the value of the slope parameter b . The corresponding values of 
our low-energy parameters are |13| : 



a) aw = 3.35 b) a nn = 4.10 , (62) 

respectively, while f3 n7T = 1.17 in both cases. The isospin 1 = 2 parameters are taken as in Ref. 
$2]. Given a matching point sq, we now have all the ingredients needed for a determination of 



the constants 7/ via (|55]) and the evaluation of the dispersive amplitudes from ( |5lD 
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In Figs. 6a and 6b we compare the dispersive evaluation of the amplitude H(s) (or 
more precisely 100 M%H(s) for < s < 4M 2 and 100M%\H(s)\ above threshold) for 77 -> 
7r°7r°, corresponding to the two parameter sets ( j60j ) and ( plj ) of 7r — 7r phases to our 0(p 5 ) 
chiral calculation. The experimental points correspond to the Crystal Ball data |]] under the 
assumption that only the S-wave contributes to the cross section. One observes that in the case 
of the phases determined from flBTf ) (Fig. 6a) the two dispersive representations (corresponding 
to So = and so = 4M 2 ), the chiral calculation, as well as the experimental data are all rather 
close to each other in all the physical low-energy region of interest. On the contrary, in the case 
of the parameters ( |6TD the dispersive amplitude always exceeds both the experimental data 
and the perturbative amplitude. This is presumably due to the fact that this set of phases is 
caracterized by a slope parameter which is too large (the scattering length being chosen to be 
ao = 0.26 in both cases). We note that the 77 — > 7r°7r° data provide only indirect information 
on the low energy behaviour of the 7r — it phases. On the other hand, the chiral representation 
allows to constrain the low energy parameters and (3^, which are related in a non-trivial 
way to the phases ||13| . One further observes that the difference between the two dispersive 
curves with matching at and 4M 2 , respectively, is most significant around s — 0. The 
amplitude is rather small in this region and corrections appear to be large. For instance the 
0(p s ) contribution coming from the vector mesons is about 50% of the 0(p 5 ) chiral amplitude 
at s = , while it represents a negligible 5% at threshold. 

Finally, we turn to the question of the pion polarizabilities. Due to our lack of knowledge 
of the 0(p 5 ) corrections A/ (see Eq. fl5B])) the charged pion polarizabilities can only be given 
to generalized 0(p 4 ) order: 



a 

«7T + 



7T 7f 



A™) + t^t^^k + 4(L 9 + L10) + 0{p 5 ) 



144tt 2V "" ' 576tt 2 



(63) 



where 



M n , n w M 2 
= 2{a nK -i) " ■- ' "" 

For the neutral pion, we find 



X K = 2K K - 1)^L + (1 _ ^ K )(i + _*.) . (64) 



a 



M n F% 



(aw - 4/3^) + r } n - 2c 



288tt 2 ^ ' 576tt 2 



(65) 



The difference with the standard 0(p 4 ) situation is that, firstly, a wn and /3 nn (as well as a n x 
and (3 n K) can depart significantly from the value 1 (the values used in the table below are those 
quoted in (p2] ) ). Secondly, we have extra contributions which are of chiral order five. In order 
to get a feeling of how significant these differences are, we have collected a few numerical values 
in Table (1) below: 
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a n o 


f3 n o 


a n + 




standard 0(p A ) 


-0.50 


0.50 


2.74 


-2.74 


standard 0(j9 6 )(ll| 


-0.35 


1.50 






generalized 0(p 5 ) 


0.44 


-0.44 


3.47 


-3.47 


dispersive(s = 4M 2 ) 


-0.76 


1.78 


3.33 


-3.05 


dispersive(s = 0) 


0.96 


0.07 


3.61 


-3.33 



Table 1: Results for the electric (a) and magnetic (/3) pion polarizabilities in 

units o/10~ 4 fm 3 . 

The dispersive results are evaluated using the phase-shift parameters (|60"D, and the values of 
a W7r , P W7r used in the 0(p 5 ) expressions (63-65) correspond to the same phases (see (62a) ). We 
proceed as follows. First, the difference in the electric and magnetic polarizations are obtained 
from 



a n o — f3 n o 



-^limiMf) 

M v s=o s 



(66) 



(/at being calculated from the representation fl5T|)) and similarly for the n + with f^ replaced 
by f C - Next, we assume that the polarization sum is saturated by the spin one meson pole 
contributions (the conventions and notations follow Ref. £[T]]): 



(67) 



with C w = 0.67, C p = 0.12, C bl = 0.53 (all in GeV~ 2 ) for the vr°, while for the charged pion 
the sum is given by the same formula with C w = 0, C p = 0.06, and C bl as before. These values 
are obtained from the experimental rates V° — > 7r°7 and V + — > 7r + 7 respectively. Note the 
strong breaking of isospin symmetry for the p which is due to p — uj mixing. 

The table shows that the difference between the standard and the generalized is 
rather minor for the charged pion which, in both cases, is dominated by L 9 + L w . Things 
are quite different for the neutral pion since, in this case, we find the electric polarizability to 
be positive. This result is essentially caused by the 0(p 5 ) constant c in formula d65|). This 
contribution is much larger than the corresponding one in the 0(p 6 ) calculation of Ref. [|ll 



because of the l/(r — 1) factor in (26) which is much larger in our case. It is also striking 
that the dispersive results with different matching points to the amplitude yield rather 
different numbers for the polarizations. Note finally that the kaon loop contribution, which one 
would a priori expect to be negligible (and which is not included in the numbers shown in the 
table) would further increase the generalized 0(p 5 ) value of a n o by 20%. This suggests that 
the neutral pion polarizabilities may be as difficult to control theoretically as they are difficult 
to measure experimentally. 
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5 Summary and concluding remarks 



i) At the one loop order of generalized xPT, the 77 — > 7r°7r° amplitude consists of three 
parts which are all separately finite: The one loop 0{p i ) and 0(p 5 ) parts, and the tree level 
0(p 5 ) contribution, the latter representing a constant shift in the amplitude H(s). Neglecting 
kaon loops and Zweig rule violating effects (described in standard xPT by the constants L 4 
and Lq), the whole amplitude can be expressed in terms of the parameter and of the single 
constant c describing the tree contribution of We have derived a low-energy theorem 

which relates the constant c to the experimental cross section for e + e~ — >hadrons and to a nn . 



ii) The constant a n7r is measurable in the low-energy n — it scattering ||13| . Its leading 
0(p 2 ) part a l ^ d is a function of the quark mass ratio r = m s /rh. Within the standard xPT, a n7T 
should remain close to u l ™ d = 1 (r ~ r 2 = 25.9). The generalized xPT admits a substancially 
larger value of a n ^, typically a n7T ~ a l ^ d < 4, corresponding to r > r\ = 6.3. For energies 
< 450 MeV, the generalized one loop cross section agrees with the Crystal Ball data []J 
within errors, provided that a nn > 3 (implying r < 10). The low-energy cross section is barely 
affected by the uncertainties in the determination of the constant c. 



iii) At order 0(p 5 ) of generalized xPT, the neutral pion polarizability (a — /3) no becomes 
positive, typically (a — j3) wo = (1.04 ± 0.60) x 10~ 4 fm 3 for a nn = 3. This has to be compared 
with the values —1 x 10~ 4 fm 3 and —1.85 x 10 _4 fm 3 predicted by the standard 0(p 4 ) and 0(p 6 ) 
orders [|ITJ, respectively, (a + (3 = up to and including order 0(p 5 ).) This result is due to 
the fact that in generalized %PT the negative 0{p A ) contribution is strongly suppressed due 
to the proximity of the Adler zero, which for a l ™ d — > 4 moves towards the point s = 0. The 
polarizability is then dominated by the positive 0(p 5 ) tree contribution. On the other hand, 
our prediction for the charged pion polarizabilities (Table 1) does not differ very much from the 
standard 0(p 4 ) result. Notice that in this case, we have no quantitative control of the 0(p 5 ) 
constant contribution. 



iv) To illustrate the convergence rate of the generalized %PT, one may first compare 
the 0(p 4 ) and 0(p 5 ) contributions to the amplitude H(s). For a nn = 3, the value of H(s) 
(multiplied by lOOM^) at threshold is 8.64, of which 7.23 come from 0(p 4 ) loops, 0.81 from 
0(p 5 ) loops and the remaining 0.60 represent the 0(p 5 ) tree contribution. Next, in order to 
control whether, at sufficiently low energies, the one loop generalized correctly accounts 
for rescattering effects and final state interactions, we have compared our 0(p 5 ) perturbative 
amplitude with a dispersion-theoretic amplitude which satisfies exact S-wave unitarity. For 
> 3 and for 7r — 7r phase-shifts not characterized by a too large slope parameter, we have 
found that the dispersive and perturbative amplitudes do indeed agree reasonably well for 
E n7T < 450 MeV, and that there is no conflict with unitarity in this energy range. These facts 
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enforce the hope that higher orders of generalized xPT will remain sufficiently small, such as 
not to spoil the agreement with experiment. 



It might be legitimate to ask that, at sufficiently low energies, xPT should correctly repro- 
duce experimental data without resorting to higher orders, unless it is justified by kinematical 
reasons or by the high precision of the data. In the case of the reaction 77 — > 7r°7r°, the standard 
expansion of C e ^' fails to satisfy this requirement. The present analysis suggests a possible 
interpretation of this fact: At each order, the standard xPT misses important symmetry break- 
ing terms which are unduly relegated to higher orders. In order to recover these terms and to 
reach agreement with experiment, one then has to go up to unnaturally high orders of xPT- 
If this conjecture turned out to be correct, similar phenomena could be expected elsewhere 
(low-energy 71 — 71 scattering, pion scalar form factors, symmetry breaking aspects of de- 
cays, etc.). In order to decide whether the generalized xPT provides a relevant improvement of 
the standard expansion of C e ^', and whether the values of the various low-energy parameters 
(-Bo, t = m s /m,...) should be revised correspondingly, additional experimental informations 
are needed. More precise data on low-energy 77 — > 7r°7r° and on the related process rj — > 7r°77 
will be welcome. Decisive information might come from the model independent determination 
of low-energy tx — 71 phase shifts and of Kn form factors at Da$ne, as well as from a direct 
measurement of light quark masses in exclusive tau decays p5 . 
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Appendix 

We have gathered, in this Appendix, a few formulae related to the discussion in Section 2. 
The explicit expressions of the combinations M 2 a Tn and (M K — M T ) 2 + 2M v M K a^K in terms 
of the low-energy constants of and of read as follows: 

MX- = 2mB + 4rh 2 (4A + 8ZS + rZ$) 

- 4m 2 5 [3£ + |(6 + r)] 

- 16m 3 A [3£ + 2£(3 + r)] 

- 8m 3 Z 5 [3£(4 + r) + £(24 + 14r + r 2 )] 

+ m 3 [81 Pl +p 2 + 2p 4 (82 + 16r + r 2 ) + p 5 (2 + r 2 ) + 12p 7 (2 + r)(14 + r)] , 
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and 



(M K - M w f + 2M n M K a wK = 

= mB {3 + r)+ m 2 A {17 + 14r + r 2 ) + 2m 2 Z 5 (18 + 17r + r 2 ) 

- m 2 £ [f (1 + r)(ll + r) + 2|(12 + 17r - r 2 )] 

- m 3 A [£(l + r)(29 + 18r + r 2 ) + 2£(40 + 71r + 18r 2 - r 3 )] 

- m 3 Z 5 [(l + r)(68 + 50r + 2r 2 ) + 2£(2 + r)(64 + 42r - 2r 2 )] 

3 

+ 2m 3 [ -pi (43 + 50r + 14r 2 + r 3 ) 
+ ^P 2 (3 + r 2 ) 

+ ip 4 (146 + 188r + 59r 2 + 3r 3 ) 
2 

+ Ip 5 (2 + r 2 )(3 + r) 

+3p 7 (2 + r)(30 + 29r + r 2 )] . 

The order 0(p 3 ) expressions for the pseudoscalar masses are given as: 

M 2 = 2mB + Am 2 A + Am 2 Z^(2 + r) 

- 4m 2 5 [£ + e(2 + r)] 

- 8m 3 A [£ + £(2 + r)] 

- 8 m 3 (2 + r)Z 5 [£ + e(2 + r)] 

+ m 3 [ 9pi + p 2 

+2p 4 (10 + 4r + 4r 2 ) + p 5 (2 + r 2 ) 
+ 12p 7 (2 + r) 2 ] , 

and 

Ml = mB (l +r) +m 2 A (l + r) 2 + 2m 2 Z 5 (2 + r)(l + r) 

- m 2 (l + r)B [£(l + r) + 2£(2 + r)] 

- m 3 (l + r) 2 A [£(l + r) + 2£(2 + r)] 

- 2m 3 (l + r)(2 + r)Z 5 [£(l + r) + 2#(2 + r)] 

3 1 
+ m 3 (l + r) [ -pi(l + r + r 2 ) + -p 2 (l - r + r 2 ) 

+3p 4 (2 + 2r + r 2 ) + ^p 5 (l + r + r 2 ) 
+6p 7 (2 + r) 2 ] . 

Combining these formulae, one may obtain the expressions for a n7T and of a w K- At leading 
order, i.e. when £, £, p lr ..p 7 are set equal to zero, one can express mfi an d m 2 A in terms of 
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the quark mass ratio r and of the pseudoscalar masses 

Ml 



m 2 A 



2(r 2 - 1) 
Ml 



[(r - n)(r + n + 2) - 2C(r 2 - r)(2 + r) 



2r 5 



[r 2 - rj 



from which one infers the expressions (f[8|) for a l ^ d and for at^jf 
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Figure captions 



Figure 1: The vertices from + O 3 ^ and from O 5 ^ (for the last one) which enter the 
computation of the order 0(p 5 ) amplitude H(s). In the penultimate and antepenultimate 
graphs, the mass-shell conditions p\ = p\ = M% for the outgoing tt° lines were used. The 
constants a nn , a^K, A™, PnK, as well as the constant c which appears in the last vertex, are 
defined in the text and in the Appendix. 

Figure 2: The cross section (7(77 — > 7r°7r°, | cos#| < Z) as a function of the center of mass 
energy E at order 0(p 4 ) for r = 25.9 (dotted line), r = 10 (dash-dotted line), and for r = 6.3 
(solid line). The data points are taken from Ref. and the vertical error bars correspond to 
the statistical errors alone. 

Figure 3: The cross section (7(77 — > tt°tt°, | cos#| < Z) as a function of the center of mass 
energy E at order 0(p 5 ) for a n7r = 1 (dotted line), a n7r = 3 (dash-dotted line) and for a nn = 4 
(solid line), with £/£ = in all three cases. 

Figure 4: The cross section 17(77 > vr 7r°, | cos6\ < Z) as a function of the center of mass 
energy E for a T7T = 3 (solid line) and the same cross section, but with c set equal to zero 
(dotted line). In both cases, £/£ = 0. 

Figure 5: The uncertainty in the order 0(p 5 ) cross section 17(77 ~~ > vr°7r°, |cos#| < Z), for 
a nn = 3, coming from the variation of the Zweig rule violating parameter £/£ between -0.2 
(lower dotted curve) and +0.2 (upper dotted curve). The solid curve corresponds to £/£ = 0. 

Figure 6: The order 0{p 5 ) amplitude 100M||#(s)|, for s > 4M%, &ndl00M*H{s), for < 
s < 4M^), shown by the full line, as compared to the dispersive results obtained by matching 
with the xPT expression at s = 4M^ (dotted curve) or at s = (dash-dotted curve) for the set 
of phase shifts given by the parameters of Eq. (|60|) (Fig. 6a) and of Eq. (|B"TD (Fig. 6b). 
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